Universal distribution of transparencies in highly conductive Nb/AlO^/Nb junctions 
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We report the observation of the universal distribution of transparencies, predicted by Schep 
and Bauer [Phys. Rev. Lett. 78, 3015 (1997)] for dirty sharp interfaces, in uniform Nb/AlO^/Nb 
junctions with high specific conductance (10 8 Ohm -1 cm -2 ). Experiments used the BCS density 
of states in superconducting niobium for transparency distribution probing. Experimental results 
for both the dc I — V curves at magnetic-field-suppressed supercurrent and the Josephson critical 
current in zero magnetic field coincide remarkably well with calculations based on the multimode 
theory of multiple Andreev reflections and the Schep-Bauer distribution. 



The basic characteristic of a mesoscopic conductor is 
its set of transmission coefficients, or "transparencies" , 
defined as the eigenvalues of the transfer matrix connect- 
ing all incoming electronic modes to outgoing modes (for 
a thorough review, see Ref. The set of transparen- 
cies for a given conductor determines all its transport 
properties including dc current and broadband current 
noise. 

In 1982 Dorokhov showed Q that the distribution of 
transparencies in diffusive conductors is universal, i.e., 
does not depend on dimension, geometry, carrier density, 
and other sample-specific properties: 
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where G is the average conductance and Go = 2e 2 //i. 
The universality of Eq. (Q) is responsible in particular 
for the universal value Sj(0) = 2e//3 of shot noise in 
diffusive conductors ||]4| (here / is the dc current). This 
suppression of shot noise in comparison with its Schottky 
value has been observed experimentally J 5j-|7[ and may 
serve as an indirect confirmation of Eq. (|l). 

Recently, Schep and Bauer showed |§) that the distri- 
bution of transparencies of a disordered interface is also 
universal, but is given by an expression different from 
Eq. (Eh g: 
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(This distribution leads to a shot noise value of Si(Q) = 
lei /I Q). In this work we report a strong experimental 
evidence that the transparency distribution in sub-nm- 
thick aluminum oxide barriers is very close to Eq. (j^) 
while being substantially different from Eq. (Q). 

The determination of the transparency distribution 
may be assisted by the fact that due to the BCS sin- 
gularity in the density of states at the edges of the su- 
perconductor energy gap 2A(T), both the Cooper-pair 
and quasiparticle transport through Josephson junctions 
are highly nonlinear. In particular, quasiparticle transfer 
strongly increases at the "gap voltage" V g = 2A(T)/e, 



while below this threshold the transport is dominated by 
multiple Andreev reflections (MAR) [|ll]-[l3| , resulting in 
a pronounced sub-harmonic gap structure at V n = V g /n. 
This structure is very sensitive to the number and trans- 
parencies of the modes, and has been successfully used 
]l4] [i"6| to determine transparencies of atomic-size point 
contacts with a few propagating modes. However, as will 
be shown below, the fact that various features of MAR 
transport (the current jump at V = V g , the excess cur- 
rent at V > V g , and subharmonic structure at V < V g ) 
are sensitive to different ranges of the transparency distri- 
bution also allows probing of the distribution in junctions 
with a much larger area (and hence a very large number 
of propagating modes). 

Our samples were made using in-situ-fabricated 
Nb/A10 x /Nb trilayers which were deposited on oxidized 
Si wafers in a cryopumped vacuum system with a base 
pressure of 5 x 10~ 8 Torr. A dc magnetron sputtered 
150-nm niobium base electrode was covered by a 8-nm- 
thick aluminum film (also using dc sputtering) . Without 
breaking the vacuum, an AlO^ layer was formed by ther- 
mal oxidation with well controlled dry oxygen exposure. 
For the junctions discussed below (specific normal con- 
ductance close to 10 8 Ohm _1 cm -2 ), the oxidation was 
carried out at 1.0 mTorr of O2 for 10 minutes. After 
the trilayer deposition was completed by sputtering of a 
150-nm-thick Nb counter electrode, junctions of various 
areas (from 0.25 to 1 /im 2 ) were formed by e-beam pat- 
for details, 
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The resulting junctions 
were highly uniform; for example, for 1 /im 2 junctions 
the full spread of normal conductance was below ±3%; 
a slightly larger (±8%) spread of deep-submicron junc- 
tions may be readily explained by their area definition 
uncertainty. The junction homogeneity was further con- 
firmed by the fact that the Josephson supercurrent could 
be almost completely suppressed by a magnetic field cor- 
responding to the insertion of one flux quantum into the 
junction - see Fig. [l]. 

The solid lines in Fig. [l] show the experimental dc I— V 
curves of two junctions of different area at a bath temper- 
ature of To = 1.8 K, with the Josephson critical current 
suppressed by a magnetic field parallel to the film plane, 
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while the solid lines in Fig. g show the measured differ- 
ential conductances of the same samples. These curves 
exhibit well pronounced sub-gap structure, indicative of 
MAR transport. However, the conduction peak posi- 
tions V n at higher voltages deviate gradually from the 
expected dependence 1/n, indicating self-heating of the 
samples. This heating was higher in larger junctions and 
may be calculated with good accuracy using a simple 
model where the main temperature gradient between the 
self-heated junction and the bath is parallel to the sub- 
strate, in the niobium electrodes adjacent to the junction; 
from this picture one should expect overheating to scale 
roughly as the junction area divided by its perimeter, in 
agreement with experiment. 
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FIG. 1. DC I — V curves for two samples of area (a) 
0.25/im 2 and (b) 1.0/im 2 . Solid lines: experimental results 
normalized with 1/G = R N = R c ^ pt ■ Dashed lines: MAR the- 
ory using Eq. (|), (a) R N = 1.01R e * pt , (b) R N = 1.04R^ pt . 
Dotted lines: MAR theory using Eq. ([[]), (a) R N = 0.82ff 
(b) R N = 0.84R^ pt . 
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The points in Fig. ^ show the dependence of the energy 
gap A(T), where T is the junction temperature includ- 
ing self-heating, on the power P = IV dissipated in the 
junction, read off from the peak positions. Assuming 
the BCS temperature dependence of the energy gap, the 



data show that at V = V g (T) the junction temperature 
is close to 6.2 K and 7.6 K for, respectively, the 0.25^m 2 
and 1.0/mi 2 junctions. The solid and dashed lines in 
Fig. H show the theoretical dependence of the gap on P 
assuming a specific heating model. In this model we ap- 
proximate the experimental temperature dependence of 
the heat conductance k(T) of superconducting niobium 
p9[ with a one-parameter parabolic expression. The re- 
sulting dependences provide a reasonably good, smooth 
interpolation of the experimental heating data. 
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FIG. 2. Differential conductance for the same samples, and 
with the same fitting values of Rn, as in Fig. 

Dashed curves in Figures |l| and || show the theoreti- 
cal dependences which were obtained by averaging the 
results of the MAR theory jl2|,[l3| over the distribution 
of transparencies given by Eq. (^) |2(if| . Apart from the 
incorporation of the heating model described above, two 
additional minor adjustments were made when calculat- 
ing these curves. First, the assumed value of normal 
resistance Rn was allowed to differ by a few percent 
(within experimental uncertainty pl|) from the exper- 
imentally determined value R c ^ pt . Second, we have in- 
troduced a Gaussian distribution of A(T) with an r.m.s. 
spread of 2%. Such a spread is typical for any super- 
conductor junction and may be readily attributed to the 
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anisotropy of the superconducting gap in polycrystalline 
electrode films. 



(middle curve) [p3fl . It is evident that the experimental 
data agree with the Schep-Bauer distribution. 
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FIG. 3. Dependence of the superconducting gap A(T) on 
the dissipated power P = IV for the same two samples as in 
Figs. [l] and |^. Points: fit to positions of the peaks in Fig. |^ 
(see the text). Lines: Calculated A(T) based on our heating 
model and the BCS gap temperature dependence. 

We believe that the agreement of the experimental 
data and theoretical curves based on the Schep-Bauer 
distribution is remarkable. In order to see that this 
agreement could not result from the fitting procedure 
described above, dotted curves in Figs. [l] and || show the 
results of our best attempt to fit the data with the MAR 
theory results averaged using the Dorokhov distribution 
(with a similar account of self-heating) . It is evident that 
these curves are rather far from the data. Moreover, in 
this best fitting attempt we have selected the values of 
Rn rather distant from R c ^ pt (see the Fig. [j] caption); if 
the latter values were used, the theoretical lines would 
pass considerably higher than the experimental plots. 
This is immediately visible from the values of the "ex- 
cess current" defined as I cxc = I(V) — GV at V ^> V g . 
Averaging the results of Rcf. fl2|| for T = with the 
Schep-Bauer distribution Eq. (H), we get 
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while for the Dorokhov distribution the excess current is 
substantially higher, I cxc w 1.467GA(0)/e ||). 

Another (though less spectacular) evidence of the va- 
lidity of the MAR theory combined with the Schep-Bauer 
distribution comes from the temperature dependence of 
the Josephson critical current (in zero magnetic field), 
which is also sensitive to the mode transparency distri- 
bution. Points in Fig. |] show the measured temperature 
dependence of the critical current in the same samples 
as in Figs. [I] ||. Lines show the Ambegaokar-Baratoff de- 
pendence for tunnel junctions (upper curve), the MAR 
theory using the Dorokhov distribution (lower curve), 
and the MAR theory using the Schep-Bauer distribution 
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FIG. 4. Temperature dependence of the critical current I c 
in the two samples of Figs. Open circles: A = 0.25 /i 2 ; 

filled circles: A — 1.0 /i 2 ; dashed-dotted line: the Am- 
begaokar-Baratoff result for tunnel junctions; dashed line: 
MAR theory using Eq. (|^); dotted line: MAR theory using 
Eq. §. 

The absolute value of the low-temperature critical 
current is different in all three models used in Fig. |J: 
/ c (0)/(GA(0)/e) = 1.57; 1.92; 2.08 for, respectively, a 
tunnel junction, disordered barrier with the Schep-Bauer 
distribution of transparencies, and SNS junctions with 
Dorokhov's distribution ^6|. The observed value lies 
within a few percent of the prediction for the Schep- 
Bauer distribution; the difference can again be attributed 
to the error in the normal resistance definition |^l| . 

The fact that transport in disordered A10 x barriers of 
finite thickness (of the order of 1 nm, i.e., much thicker 
than the Fermi wavelength Xp in the junction electrodes) 
is described by Eq. (^|) so well may seem rather surpris- 
ing, since its derivation in Ref. || assumed that the bar- 
rier is a strongly-disordered region with thickness d much 
smaller than Xp. However, this distribution may be de- 
rived from a different model which does not rely on this 
assumption. 

It is well known that resonant tunneling through a sin- 
gle localized site leads to the following transparency: 
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where e is the state energy, ep is the Fermi level in the 
junction electrodes, T is the tunneling width for a site 
in the barrier center, x is the site deviation from the 
center, and a is the localization radius. If T is so large 
that the first term in the denominator is unimportant, 
for a system with a uniform spatial distribution of sites, 
Eq. (Q) immediately gives the Dorokhov distribution. On 
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the other hand, if x ~ while the spread of e is much 
broader than T, Eq. (0) yields the distribution (||). This 
condition applies to strongly disordered barriers like ours 
p7| if their thickness d is smaller than a (though possibly 
much larger than A^) and their spread of atomic local- 
ization energies is larger than T. 

In summary, the excellent agreement between the ex- 
perimental data and the results of the MAR theory com- 
bined with the Schep-Bauer distribution of transparen- 
cies provides very strong evidence that our Nb/AlO^/Nb 
junctions are well described by this distribution, while 
being far from, e.g., the Dorokhov distribution. Since the 
distribution (||) has the same universal nature as Eq. (Q), 
this result is of considerable general importance for meso- 
scopic physics [p8| . 

Our result is also of substantial importance for applica- 
tions. It shows that transport in niobium-trilayer Joseph- 
son junctions with high specific conductance and hence 
high critical current density (up to at least 200 kA/cm 2 ) 
may be due to a fundamental mechanism rather than rare 
defects such as pinholes, etc. This gives every hope that 
these junctions may be even more reproducible than in 
our first experiments: assuming that the barrier trans- 
parency is only correlated at distances of the order of 
its thickness (about 1 nm), we may estimate that the 
minimum r.m.s. spread of the critical current is below 
1% even for deep-submicron junctions. Together with 
the fact that such junctions are intrinsically overshunted 
, this makes them uniquely suitable for several impor- 
tant applications in superconductor electronics, including 
ultrafast digital RSFQ circuits of very high integration 
scale - see, e.g., Ref. pjj . 
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